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Abstract.  Cementitious  materials  have  material,  geometric  features  ranging  from 
molecular/nano,  micro-,  tneso-  and  macro  scales.  Features  and  changes  in  material 
chemistry/nano  scale  influence  the  hydration  process,  formed  micro  scale  morphology, 
associated  properties  and  behavior  at  engineering  length  scales.  In  the  present  paper,  effect  of 
material  ion  exchange  on  mechanical  stiffness  properties  of  nano  level  hydrated  cement  paste 
constituent  -  calcium  silicate  hydrate  (C-S-H)  based  on  its  material  chemistry  structure  are 
studied  following  a  molecular  dynamics  (MD)  computational  modeling  methodology. 
Calcium  ions  are  replaced  with  Magnesium  ions  in  traditional  C-S-H  Jennite  crystal  structure 
of  the  hydrated  cement  constituent.  Traditional  C-S-H  Jennite  structure,  and  Magnesium 
modified  structures  obtained  by  calcium  exchange  of  magnesium  are  used  in  the  MD  based 
material  modeling  analysis.  Extensive  modeling  analysis  established  optimal  computational 
MD  modeling  parameters  and  molecular  material  sizes  employed.  The  present  study  clearly 
shows  the  influence  of  material  chemistry  due  to  Magnesium  ion  exchange  in  CSH-Jennite 
with  predicted  elastic  modulus  that  changes  with  the  level  of  Magnesium  in  CSH-Jennite 
material  chemistry  structure. 

1  INTRODUCTION 

Cementitious  materials  have  material,  geometric  features  ranging  from  molecular/nano, 
micro-,  meso-  and  macro  scales.  Features  and  changes  in  material  chemistry/nano  scale 
influence  the  hydration  process,  formed  micro  scale  morphology,  associated  properties  and 
behavior  at  engineering  length  scales.  Cement  paste  is  made  by  combining  dry  Portland 
cement  and  water.  The  process  of  making  cement  paste  is  referred  to  as  hydration,  during 
which,  there  is  loss  of  workability,  solidification  and  hardening.  The  hydration  of  the  different 
clinker  phases  results  in  different  hydration  product  phases,  which  makes  up  the  cement  paste. 
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The  hydrated  cement  phases  include;  calcium  silicate  hydrate  (C-S-H),  calcium  aluminum 
hydrate  (C-A-H)  and  calcium  hydroxide  (C-H)  '1>2'.  In  cement  paste,  four  levels  of  scales  have 
been  identified13'  i.e.  molecular  /nano-scale  level  determined  by  the  material  chemistry 
features,  micro-scale  with  microstructural  morphology  features,  meso-scale  and  macro-scale. 
Distinct  length  scales  and  features  are  based  on  the  requirement  that  each  scale  should  be 
separated  from  the  next  scale  by  at  least  one  order  of  length  magnitude,  which  is  also  stated  in 
the  literature  as  a  prerequisite  for  the  application  of  continuum  micromechanics'-1’  2'.  As  a 
result  of  the  multi-scale  nature  of  cement  paste,  the  properties  and  material  composition  of  the 
macro  and  engineering  scale  cement  can  be  controlled  or  affected  by  the  molecular  level.  A 
potential  ability  to  alter  and  manipulate  the  chemistry  level  structure  of  hydrated  component 
of  traditional  C-S-H,  a  structure  predominantly  influenced  by  calcium  -  a  major  constituent  of 
the  cement  paste  at  the  fundamental  material  chemistry  level  -  will  imply  that  the  elastic  and 
deformation  properties  of  the  paste  would  be  impacted.  Figure  1  shows  the  different  levels  of 
scale  associated  with  cementitious  materials. 
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Figure  1:  Schematic  of  multi-scale  structure  of  cement  based  materials  [3] 


Traditional  cement  paste  is  primarily  calcium-based,  but  there  is  interest  in  the  engineering 
community  for  alternate  formulations  of  cement;  which  may  provide  tailored  properties  for 
specific  applications,  could  utilize  local,  and  cheaper  materials.  Another  factor  influencing  the 
alternate  forms  cement  cited  in  the  literature  is  the  need  to  reduce  the  energy  required  to 
produce  the  clinker  phase,  and  the  carbon  foot  print  of  the  cement  industry"1'.  Hence  there 
have  been  efforts  to  look  for  alternate  cement  formulations  that  are  not  calcium-based,  might 
be  cheaper,  and  have  the  potential  for  reducing  carbon  footprint  in  the  manufacturing  stage  of 
calcium  based  cement.  Several  forms  of  cement  formulations  has  been  sought,  for  instance, 
other  fonns  of  belitic  cements  are  sought  by  trial  and  error  to  improve  the  reactivity  of  the 
clinker  phase.  Some  additional  approaches  that  have  been  used  trying  to  improve  the 
properties  of  cement-based  materials  include  thermal  processing'4’5'  and  addition  of  new 
chemical  compounds  and  ions  '4’5’6'.  Most  common  types  of  ionic  substitution  which  has  been 
used  in  cement  include;  Mg2+  for  Ca2+,  2A13+  and  2Fe3+  for  3Ca2+,  and  2Ca2+  for  Si4+  [6’7].  The 
present  work  investigates  the  influence  of  replacement  of  calcium  ion  with  magnesium  ion  in 
the  material  chemistry  level  of  C-S-H  with  a  focus  on  mechanical  stiffness  detennined  via 
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computational  material  modeling. 

Several  studies  have  focused  on  different  aspects  of  cement  material  improvement.  Hegoi 
and  co-workers*6'  studied  the  effects  of  the  presence  of  chemical  substitution  on  the 
physicochemical  properties  of  cement  clinker  phases  -  alite  and  belite.  They  incorporated 
Mg2+,  Al3+  and  Fe3+  into  the  structure  using  classical  forcefield  methods  and  reported  that  the 
crystallographic  site  within  the  unit  cell  is  equally  probable  for  Mg,  A1  and  Fe  substitution.  It 
has  also  been  reported  that  Mg  incorporation  does  not  change  the  electronic  structure 
appreciably  compared  to  A1  and  Fe  incorporation*6',  and  that  there  are  no  preferential 
substitution  for  any  calcium  sites,  when  Mg2+,  Al3+  and  Fe3+  were  incorporated.  The  main 
focus  of  their  research  was  on  the  reactivity  of  the  clinker  phase  for  reduced  energy 
production. 

As  discussed  in  the  literature,  Mg  provides  a  feasible  incorporation  option  in  cement  paste. 
In  a  study  of  the  hydration  process  of  cement,  Stephan  and  co-workers*7-',  provided  analytical 
evidences  that  in  low  concentrations,  Mg"  does  not  change  the  hydration  process  of  calcium 
trisilicate  (C3S)  -  a  component  of  cement,  although  the  reactivity  may  increase  with  weight 
percent  and  with  cement  age.  Also,  because  of  the  relative  availability  and  the  ionic 
similarities  to  Calcium,  Mg  can  be  considered  as  a  potential  candidate  for  ionic  replacement. 
Additionally,  with  the  replacement  of  Ca  by  Mg,  a  lower  temperature  for  processing  the 
clinker  phase  can  be  achieved. 

The  focus  of  this  paper  is  to  understand  the  influence  of  the  effects  of  the  exchange  of 
Calcium  ion  with  Magnesium  ions  in  hydrated  cement  material  chemistry  structure  based  on 
the  mechanical  stiffness  properties  of  the  modified  material  chemistry  configuration  using 
computational  material  modeling  and  simulations.  Different  computational  simulation 
methods  have  been  employed  in  the  past  to  study  heterogeneous  materials.  Manzano  et  al'6' 
used  a  combination  of  forcefield  and  DFT  atomistic  simulations,  while  others  have  used  ab 
initio  methods  to  study  the  structure  of  C-S-H*8,9'.  There  are  also  prior  studies  on  the 
mechanical  behavior  of  heterogeneous  materials  using  multi-scale  modeling,  and  studies  to 
understand  cement  hydration  with  finite  element  modeling,  as  well  as  temperature 
dependency  of  the  microstructure  of  cement  hydrates*3’10,11'  reported  in  the  literature.  The 
present  paper  presents  and  discusses  results  from  our  work  on  material  chemistry  modeling  to 
understand  the  effect  of  Calcium  ion  substitution  with  Magnesium  on  chemistry  level 
structure  of  C-S-H  on  the  predicted  mechanical  stiffness. 

2  MATERIAL  CHEMISTRY  LEVEL  MODELING  OF  C-S-H 

Material  chemistry  level  modeling  following  the  principles  and  techniques  commonly 
grouped  under  Computational  Material  Science  is  one  of  the  key  facilitator  of  material  science 
that  is  experiencing  fast  growth  pace  and  forms  the  basis  of  the  present  work*-13'.  The  major 
hydrated  component  of  cement  paste  C-S-H  gel  is  responsible  for  the  load  bearing  attributes 
and  strength  of  cement,  and  has  a  complex  material  chemistry*-14'.  The  molecular  structural 
representations  of  C-S-H  gel  are  still  inconclusive  with  several  structural  representations  cited 
in  literature*15*.  These  include: 

•  Wollastonite  group  which  comprises  of  Foshagite  (CaXSfiOgXOHX)*16',  Hillebrandite 

(Ca2(Si03)(0H)2)[17],  Xonotlite  (Ca6Si6Oi7(OH)2)[18],  Okenite 

([Cag(Si60i6)(Si60i5)2(H20)6]4  [  Ca2(H20)9)'3H20])*19',  and  others  that  are  compiled  by 
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Richardson1151. 

•  Tobennorite  group  which  comprises  of  Clinotobennoritec  (CasSieOn  5H20), 
Clinotobermorited  (Ca5Si60i7  5H20),  Clinotobennorite  9AC  (Ca5Si60i6'(0H)2), 
Clinotobermorite  9  Ad  (Ca5Si6Oi6'(OH)2)[20’21],  anomalous  and  nonnal  Tobennorite  11  A0 
(Ca4Si60i5(0H)2'5H90  and  Ca45Si60i6'(0H)2'5H20  respectively)[22‘23]  and  Tobennorite 
Uk™. 

•  Jennite  group  comprising  Jennite  and  Metajennite  (CaySieOixlOH^  SfhO).  A  complete 
composition  and  crystal  information  for  most  of  the  C-S-H  and  other  related  phases  have 
been  tabulated  by  Richardson1 1 A 

From  various  material  chemistry  molecular  structures  of  C-S-H,  Tobennorite  14A  and 
Jennite  structures  are  traditionally  accepted,  widely  used,  and  adapted  material  molecular 
structural  representation  of  C-S-H  cement  pastef2’14’15’25'261.  In  the  present  paper,  the  structural 
configuration  of  C-S-H  Jennite  is  employed  to  study  and  understand  the  influence  of  material 
chemistry  changes  due  to  Magnesium  replacement. 

Jennite  is  a  representative  mineral  fonn  of  calcium  silicate  hydrate  (C-S-H).  It  has  a 
chemical  formula  of  Ca9Si6018(0H)6. 8(H20)[15,24,25].  The  molecular/chemistry  level 
structure  of  C-S-H  Jennite  is  defined  by  a  triclinic  unit  cell  of  dimensions  a  =  10.6  A,  b  = 
7.3  A,c  =  10.9  A,  and  angles  a  —  101.3  °  ,  /3  —  97.0  ,  y  =  109.7°.  In  Figure  2  a  crystal  of 
C-S-H  Jennite,  consisting  of  64  unit  cells  in  a  4x4x4  arrangement,  is  illustrated.  The 
molecular  structure  of  C-S-H  Jennite  is  layered,  fonned  by  sheets  of  calcium  oxide  connected 
to  short  silica  chains,  three  silica  monomers1  long,  that  are  not  connected  to  each  other.  The 
calcium  oxide  sheets  are  linked  to  each  other  by  additional  calcium  octahedral  sites.  The 
layers  containing  the  calcium  octahedral  sites  connecting  the  calcium  oxide  sheets  also 
contain  most  of  the  free  water  molecules  in  the  structure'-25  . 

The  atomic  positions  of  the  initial  unit  cell  of  C-S-H  Jennite  were  obtained  from  American 
Mineralogist  Crystal  Structure  Database  (AMCSD),  which  corresponds  to  the  C-S-H  Jennite 
structure  defined  by  Bonaccorsi  et  al}25\  The  molecular  dynamics  modeling  analysis  was 
perfonned  using  the  molecular  dynamics  analysis  code  Accelrys  Material  Studio®,  employing 
the  Discover  module  with  the  associated  atom  interaction  energy  represented  by  COMPASS 
force  field'28l 

Material  chemistry  molecular  models  of  C-S-H  are  employed  based  on  a  Molecular 
Dynamics  (MD)  modeling  methodology  to  understand  the  influence  of  Magnesium  ion 
exchange  are  based  on  a  64  unit  cell  configuration  shown  in  Figure  2.  The  material  molecular 
chemistry  system  size  of  64  unit  cells  and  other  MD  modeling  parameters  were  detennined 
after  a  careful  study'29".  Other  MD  simulation  analysis  parameters  employed  were  a  dynamic 
analysis  time  of  200  ps;  periodic  boundary  conditions  applied  to  the  molecular  system, 
temperature  and  pressure  control  based  on  Nose  Thennostat  and  Parinello  Barostat,  with  a 
time  integration  time  step  of  1  fs  (femto  second)  following  a  Velocity- Verlet  method. 


1  A  silica  monomer  is  a  single  SiCB  unit. 
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Figure  2:  Molecular  structure  of  nanoscale  C-S-H  Jennite.  Both  full  3D  atomic  structure  and  configuration 
with  Calcium  Oxide  molecules  represented  in  green  and  Si02  molecules  represented  in  yellow  are  shown 

2.1  Magnesium  Ion  Exchange  in  C-S-H  Jennite 

The  relevance  and  applicability  of  Magnesium  ion  exchange  in  hydrated  cement  C-S-H 
structure  was  discussed  earlier.  Following  this,  Magnesium  (Mg)  ion  exchange  in  traditional 
C-S-H  Jennite  was  obtained  by  replacing  the  Calcium  (Ca)  ions  in  the  molecular  structure  of 
Jennite  C-S-H  structure  with  magnesium  ions.  The  unit  cell  configuration  of  C-S-H  Jennite 
consists  of  9  calcium  ions.  The  magnesium  modified  C-S-H  Jennite  structure  based  on  a  same 
crystal  configuration  is  obtained  by  replacing  one  or  more  calcium  atoms  by  magnesium 
atoms.  The  different  magnesium  modified  C-S-H  Jennite  structures,  and  their  corresponding 
magnesium  percentage  based  on  calcium  and  magnesium  molar  mass  considered  in  the 
present  work  are  tabulated  in  Table  1|291. 


Table  1:  Number  and  percentage  of  Mg  in  ion  exchange  in  one  unit  cell  C-S-H  Jennite  configuration 


Structure 

Ca  atoms 

Mg  atoms 

Magnesium  percentage  by 
weight 

JO 

9 

0 

0.0 

J1 

8 

1 

7.0 

J2 

7 

2 

14.8 

J3 

6 

3 

23.0 

J4 

5 

4 

32.7 

J5 

4 

5 

43.0 

J6 

3 

6 

54.8 

J7 

2 

7 

68.0 

J8 

1 

8 

82.9 

J9 

0 

9 

100.0 

Different  unit  cell  configurations  of  traditional  and  magnesium  modified  C-S-H  Jennite  as 
noted  above  was  replicated  to  create  a  64  unit  cell  configuration  in  a  4x4x4  arrangement. 
Figure  3  shows  the  material  chemistry  molecular  structures  of  64  unit  cell  C-S-H  structures 
created  and  analyzed  in  the  present  work.  MD  dynamic  analysis  was  performed  to  obtain 
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associated  Virial  stresses,  potential  energy  of  the  system  and  elastic  stiffness  matrix.  The  time 
average  values  of  the  molecular  dynamic  analysis  after  reaching  a  stable  energy  configuration 
were  employed  in  detennining  these  physical  parameters. 

2.2  Predicted  Mechanical  Stiffness  Properties  from  Material  Chemistry  Structure 

Mechanical  stiffness  properties  (Bulk,  Shear,  and  Young’s  Modulus)  were  determined 
from  the  elastic  constants  obtained  from  the  second  derivative  of  potential  energy  of  the 
material  chemistry  molecular  system.  This  is  defined  by  equation  1 . 


J3:  Modified  Jennite 


Jl:  Modified  Jennite 


J2:  Modified  Jennite 


Figure  3:  Material  chemistry  structures  (initial  on  left;  final  after  MD  analysis)  of  traditional  and  Mg-modified 

C-S-H  Jennite 


1  d2U 
V  deidej 


(1) 


Here,  s,,  Sj  are  lattice  strain  components,  U  is  potential  energy,  V  is  simulation  cell  volume. 
The  6><6  elastic  stiffness  matrix,  C  can  be  obtained  from  the  above  relation. 


'C 


n 


C 


16 


Cij  = 


(2) 


vC, 


61 


c, 


66' 


Elastic  compliance  matrix,  S  =  C  can  also  be  obtained. 

(•^11  "■  *^16\ 

*  :  J 

J61  ->66' 

Hence  the  mechanical  stiffness  properties  for  the  homogeneous  stiffness  modulus  can  be 
predicted  following  Voight,  Reuss  or  Hill  average  elastic  stiffness  relations  from  the  elastic 
and  compliance  matrix  coefficients  as 


Kr 

Kv 


i 


(Sn+S22+S33)  +  2(Si2+S23+S3i) 


(Cii+C22+C33)  +  2(Ci2+C23+C3i) 
9 


(4) 

(5) 
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Kh  = 


(Kv+Kr) 

2 


Gr 


15 

4(-?ll+^22+-533)_4(-?12+^23+-53l)  +  3(-?44+^S5+-566) 


G 


V  ~~ 


(Cll  +^22  +£33  )~  (Cl2  +Q3  +C'3l)  +  3(t-44+t-’sS+^66) 
15 


Gh 


( GV+Gr ) 
2 


E  =  2G(l  +  v) 


(6) 

(7) 

(8) 
(9) 

(10) 


where;  KR=  Reuss  bulk  modulus,  Kv  =  Voight  bulk  modulus,  KH  =  Hill  average  bulk 
modulus,  Gr=  Reuss  shear  modulus,  Gv  =  Voight  shear  modulus,  GH  =  Hill  average  shear 
modulus,  E=  Young’s  Modulus,  and  C;/  and  Sy  are  the  components  of  the  stiffness  and 
compliance  matrix  respectively  ^30'31^. 


3  RESULTS  AND  DISCUSSIONS 

The  MD  analysis  methodology  and  predictive  mechanical  stiffness  values  from  the 
potential  energy  as  defined  in  section  2  were  obtained  for  traditional  and  different  Mg- 
modified  C-S-H  material  chemistry  configurations  listed  in  Table  1.  Figure  4  shows  the 
change  in  the  total  energy  for  traditional  and  Mg-modified  C-S-H  Jennite  configurations.  The 
variation  of  normalized  total  energy  is  shown  in  Figure  4. 

The  predicted  elastic  modulus  (Hill  average)  values  for  traditional  and  Mg-modified  C-S-H 
Jennite  material  chemistry  structures  are  shown  in  Figure  5.  A  sharp  increase  in  the  predicted 
elastic  modulus  is  observed  for  the  Mg-modified  C-S-H  Jennite  structure  at  33%  Mg  structure 
followed  by  a  reduction  in  predicted  elastic  modulus  values  as  the  percentage  of  Mg 
modification  further  increases.  Clearly,  these  results  show  that  predicted  elastic  modulus  have 
an  impact  with  variations  as  shown  in  Figure  5.  A  Magnesium  ion  percentage  of  33%  is 
noticed  to  provide  the  highest  increase  in  the  predicted  elastic  modulus  in  Mg-Ion  modified 
C-S-H  compared  to  traditional  C-S-H  material  chemistry  structure. 


Figure  4:  Change  in  total  energy  per  mole  for  traditional  and  Mg-modified  C-S-H  Jennite 
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Figure  5:  Predicted  elastic  modulus  (Hill  average)  for  traditional  and  Mg-modified  C-S-H  Jennite 


The  corresponding  variation  in  the  predicted  bulk  and  shear  modulus  for  various  Mg- 
modified  C-S-H  Jennite  material  chemistry  structures  is  shown  in  Figure  6.  Again  in  both 
cases,  present  results  indicate  that  there  is  a  positive  impact  with  increased  bulk  and  shear 
modulus  values  until  a  Magnesium  ion  exchange  of  about  33%. 


Mg(%) 


Figure  6:  Predicted  shear  and  bulk  modulus  (Hill  average)  for  traditional  and  Mg-modified  C-S-H  Jennite 


4  CONCLUSIONS 

Material  chemistry  and  modifications  influence  the  evolving  behavior  of  complex, 
hierarchical  materials  such  as  cement  where  material  chemistry  changes  influence  the 
hydration  process,  altering  the  microstructure  and  finally  the  engineering  scale  properties.  In 
this  paper,  mechanical  stiffness  properties  based  on  elastic  modulus  predictively  obtained 
from  traditional  and  Mg-modified  showed  that  material  chemistry  does  influence  structural 
material  stiffness  modulus  properties,  and  varies  as  the  amount  of  Mg  inclusion  increases. 
These  predicted  material  properties  are  based  on  molecular  dynamics  based  material  modeling 
with  material  chemistry  structures;  and  thus  provides  an  effective  computational  methodology 
to  understand  the  expected  variations  due  to  material  chemistry  changes  such  as  due  to 
Magnesium  ion  exchange  of  C-S-H  Jennite  discussed  in  the  present  paper.  Results  from  our 
current  analysis  is  presented  and  further  investigations  are  needed  to  explain  some  of  the  non¬ 
physical  computational  model  values  noted  in  the  predicted  values  of  bulk  modulus  that  are 


8 


Ram  V.  Mohan,  Wayne  D.  Hodo  and  Babatunde  Adebiyi 


based  on  mechanics  based  approximations  following  the  compliance  matrix  definitions  in  MD 
simulation  analysis. 

Engineering  scale  properties  of  interest  are  the  stress-strain  deformation  behavior  and 
effective  prediction  of  estimated  stress-strain  behavior  due  to  material  chemistry  changes, 
such  as  those  discussed  due  to  magnesium  exchange  in  C-S-H  Jennite  are  also  of  interest.  Our 
work  and  results  from  such  material  chemistry  level  modeling  will  be  highlighted  in  the 
presentation  and  will  be  discussed  in  future  publications. 


REFERENCES 

[1]  Van  Hoss,  H.G.  Background  facts  and  issues  concerning  cement  and  cement  data.  (2005). 

[2]  Taylor,  H.F.  Cement  chemistry.  Thomas  Telford,  (1997). 

[3]  Bernard,  O.  Ulm,  F.J.  and  Lemarchand,  E.  A  multiscale  micromechanics-hydration  model 

for  early  age  elastic  properties  of  cement  based  materials.  Cement  and  Concrete 
Research  (2008)  33(9):  1293-1309. 

[4]  Martin-Sedeno,  M.C.  et.al.  Aluminum  rich  belite  sulfoaluminate  cements:  clinkering  and 

early  age  hydration.  Cement  and  Concrete  Research.  (2010)  40(3):359-369. 

[5]  Hewlett,  P.  Lea ’s  chemistry  of  cement  and  concrete.  Butterworth-Heinemam,  (2003). 

[6]  Manzano,  H.  et  al.  Impact  of  chemical  impurities  on  the  crystalline  cement  clinker  phases 

determined  by  atomistic  simulations.  Crystal  Growth  and  Design.  (2011)  11(7):2964- 
2972. 

[7]  Stephan,  D.  Wistuba,  D.  Crystal  structure  refinement  and  hydration  behavior  of  3CaO 

Si02  solid  solutions  with  MgO,  AFO,  and  FciO,.  Journal  of  European  Ceramic  Society. 
(2006)  26(1). 141-148. 

[8]  Jansang,  B.  Nonat.  A.  and  Skibested,  J.  Modeling  of  guest-ion  incorporation  in  the 

anhydrous  calcium  silicates  phases  of  Portland  cement  by  periodic  density  functional 
theory  calculations.  In  CONMOD-2010.  3rd  RILEM  conference  on  concrete  modeling. 
(2010). 

[9]  Churakov,  S.V.  Hydrogen  bond  connectivity  in  jennite  from  abinitio  simulations.  Cement 

and  Concrete  Research.  (2008)  38(12).  1359-1364. 

[10]  Bernard,  F.  Kamali,  S.  and  Prince,  W.  3D  multi-scale  modeling  of  mechanical  behavior 
of  sound  and  leached  mortar.  Cement  and  Concrete  Research.  (2008)  38(4):449-458. 

[11]  Nakari,  K.  et  al.  Enhanced  thermodynamic  analysis  coupled  with  temperature  dependent 
microstructures  of  cement  hydrates.  Cement  and  Concrete  Research.  (2007)  37(2):  139- 
150. 

[12]  Grondin,  F.  et  al.  Multi-scale  modeling  for  the  behavior  of  damaged  concrete.  Cement 
and  Concrete  Research.  (2007)  37(  10):  1453- 1462. 

[13]  Landman,  U.  Materials  by  numbers:  computations  as  tools  of  discovery.  Proceedings  of 
National  Academy  of  Sciences  of  the  United  States  of  America.  (2005)  102(1 9):6671- 
6678. 

[14]  Nonat,  A.  The  structure  and  stoichiometry  of  CSH.  Cement  and  Concrete  Research. 
(2004)  34(9):1521-1528. 

[15]  Richardson,  I.  The  calcium  silicate  hydrates.  Cement  and  Concrete  Research.  (2008) 
38(2):  137-158. 


9 


Ram  V.  Mohan,  Wayne  D.  Hodo  and  Babatunde  Adebiyi 


[16]  Gard,  J.  and  Taylor,  H.  The  crystal  structure  of  foshagite.  Acta  Crystallographica. 
(1960)  13(10):785-793. 

[17]  Posr,  J.E.  Crystal  structure  of  hillebrandite:  a  natural  analogous  of  calcium  silicate 
hydrate  (CSH)  phases  in  Portland  cement.  American  Mineralogist.  (1995)  80:841-844. 

[18]  Hejny,  C.  and  Armbruster,  T.  Polytypism  in  xonotlite  C^SkO^  (OfTh. 
Kristallographie/International  Journal  for  Structural,  Physical,  and  Chemical  Aspects  of 
Crystalline  Materials .  (2001)  216:396-408. 

[19]  Gard,  J.  and  Taylor,  H.  Okenite  and  nekoite  (a  new  mineral).  Mineral  Mag.  (1956).  31:5- 
20. 

[20]  Merlino,  S.  Bonaccorsi,  E.  Armbruster,  T.  The  real  structures  of  clinnotobemorite  and 
tobermorite  9A  OD  character,  prototype,  and  structural  relationships.  European  Journal 
of  Mineralogy.  (2000)  12(2):41 1-429. 

[21]  Henmi,  C.  and  Kusachi,  I.  Monoclinic  tobermorite  from  fuka,  bitchu-cho,  Okoyama 
Perfecture.  Japan  J.  Min.  Petr.  Econ.  Geol.  (1989)84:374-379. 

[22]  Merlino,  S.  Bonaccorsi,  E.  and  Armbruster,  T.  Tobermorites:  their  real  structure  and 
order-disorder  (OD)  character.  American  Mineralogist.  (1999)  84:1613-1621. 

[23]  Merlino,  S.  Bonaccorsi,  E.  and  Armbruster,  T.  The  real  structure  of  tobermorite  1 1 A 
normal  and  abnormal  forms,  OD  character  and  polytypic  modifications.  European 
Journal  of  Minerology.  (2001)  13(3):577-590. 

[24]  Bonaccorsi,  E.  Merlino.  and  Kampf,  A.R.  The  crystal  structure  of  tobermorite  14A 
(plobmbierite),  a  C-S-H  phase.  Journal  of  American  Ceramic  Society.  (2005)  88(3):505- 
512. 

[25]  Bonaccorsi,  E.  Merlino,  S.  and  Taylor,  H.  The  crystal  structure  of  jennite, 
Ca9Si60i8(0H)6  8H20.  Cement  and  Concrete  Research.  (2004)  34(9):  148 1-1488. 

[26]  Gard,  J.  et  al.  A  reexamination  of  jennite.  American  Mineralogist.  (1977)  62(3-4):365- 
368. 

[27]  MinDat.org.  Jennite.  (2014).  available  from  http://www.mindat.org/min2087.html 

[28]  Sun,  H.  COMPASS:  An  ab  initio  force-field  optimized  for  condensed-phase  applications 
-  Overview  with  details  on  alkane  and  benzene  compounds.  Journal  of  Physical 
Chemistry  B.(  1998)  102(38):7338-7364. 

[29]  Adebiyi,  B.  Effect  of  material  ion  exchanges  on  the  mechanical  stiffness  properties  and 
shear  deformation  of  hydrated  cement  material  chemistry  structure  C-S-H  Jennite  -  a 
computational  modeling  study.  MS  Thesis.  North  Carolina  A&T  State  University.  (2014). 

[30]  Hill,  R.  The  elastic  behavior  of  crystalline  aggregate.  Proceedings  of  the  Physical 
Society.  Section  A.  (1952)  65(5):349. 

[31]  Liu,  Y.  et  al.  First  principles  study  of  the  stability  and  mechanical  properties  of  MC 
(M=Ti,  V,  Zr,  Nb,  Hf  and  Ta)  compounds.  Journal  of  Alloys  and  Compounds.  (2014) 
582:500-504. 


10 


